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■ 

Context. Spitzer Space Telescope observations revealed powerful mid-infrared (mid-IR) IT rotational line emission from the Stephan's 
Quintet (SQ) X-ray emitting large scale shock (~ 15 x 35 kpc 2 ) associated with a collision between two galaxies. Because H 2 forms 
on dust grains, the presence of IT is physically linked to the survival of dust, and we expect some dust emission to come from the 
\ molecular gas. 

Aims. To test this interpretation, IR observations and dust modeling are used to identify and characterize the thermal dust emission 
from the shocked molecular gas. 

Methods. The spatial distribution of the IR emission allows us to isolate the faint PAH and dust continuum emission associated with 
| the molecular gas in the SQ shock. We model the spectral energy distribution (SED) of this emission, and fit it to Spitzer observations. 

The radiation field is determined with GALEX UV, HST V-band, and ground-based near-IR observations. We consider two limiting 
i-C ■ cases for the structure of the H 2 gas. It is either diffuse, penetrated by UV radiation, or fragmented into clouds optically thick to UV. 

Results. Faint PAH and dust continuum emission are detected in the SQ shock, outside star-forming regions. The 12/24 pm flux ratio 
in the shock is remarkably close to that of the diffuse Galactic interstellar medium, leading to a Galactic PAH/VSG abundance ratio. 
However, the properties of the PAH emission spectrum in the shock differ from that of the Galaxy, which may suggest an enhanced 
+jf ■ fraction of large and neutrals PAHs. In both models (diffuse or clumpy H 2 gas), the IR SED is consistent with the expected emission 

from dust associated with the warm (> 150 K) H 2 gas, heated by a UV radiation field of intensity comparable to that of the solar 
neighborhood. This is in agreement with GALEX UV observations that show that the intensity of the radiation field in the shock is 
G uv = 1.4 ± 0.2 [Habing units]. 

Conclusions. The presence of PAHs and dust grains in the high-speed (~ 1000km s _1 ) galaxy collision suggests that dust survives. 
We propose that the dust that survived destruction was in pre-shock gas at densites larger than a few 0.1 cm 3 , which was not shocked 
£ — ' at velocities larger than ~ 200 km s . Our model assumes a Galactic dust-to-gas mass ratio and size distribution, and present data do 

£ — , not allow us to identify any significant deviations of the abundances and size distribution of dust grains from that of the Galaxy. Our 

\Q • model calculations show that far-IR Herschel observations will help constraining the structure of the molecular gas, and the dust size 

f*"^ ' distribution, and thereby to look for signatures of dust processing in the SQ shock. 

Key words. ISM: general - ISM: dust, extinction - ISM: molecules - ISM: structure - Atomic processes - Molecular processes - 
Shock waves - Infrared: ISM - Galaxies: clusters: individual: Stephan's Quintet - Galaxies: evolution - Galaxies: interactions 
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1. Introduction shocked-gas excitation diagnostics from optical emission lines 

(Xuetal., 2003). This extended region is denoted "ridge" or 

Stephan's Quintet (Hickson Compact Group HCG92, Arp 319, simply "SQ shock" in this paper. 

^ i hereafter SQ) is an extensively studied compact group of four Observations with the infrared spectrograph (IRS, 

interacting galaxies that have a complex dynamical history (e.g. H oucketal., 2004) onboard the Spitzer Space Telescope 

Allen & Sullivan, 1980; Moles et al., 1997). A remarkable fea- haye revealed a powerful mi d-infrared (mid-IR) H 2 rotational 

ture of SQ is that a giant (« 15 x 35 kpc) shock is ere- line emission from warm ^ 1Q 2 _ 1Q 3 K) mo i ecu lar gas in 

ated by an intruding galaxy, NGC 7318b (Sbc pec), collid- the SQ shock (A ppleton et al., 2006; Cluver et al., 2010). The 

ing into NGC 7319's tidal tail, at a relative velocity of ~ Hi emission is extended not only over the whole ridge, but 

1000 km s- 1 . Evidence for a group-wide shock comes from ob- in seyeral Qther structureSi inc i udmg an extension towards 

servations of an extended X-ray ridge containing shock-heated the Seyfert galaxy NGC m9> and the inter g a i actic sta rburst 

(~ 5 x 10 6 K) gas (Pietsch et al., 1997; Trinchien et al., 2003, SQ . A (Xu et aL> m9) xhe latter stracture> beyond the north- 

2005; O'Sulhvanetal., 2009), strong radio synchrotron emis- em tip of the ridgej has been shown tQ contain significant 

sion from the radio emitting plasma (Allen & Hartsuiker, 1972; CO -emitting gas (Gao & Xu, 2000; Smith & Struck, 2001; 

Sulentic et al., 2001; Williams et al., 2002; Xu et al., 2003) and Lise nfeld et al., 2002). To explain the H 2 emission from the 

SQ ridge, Guillard et al. (2009) considered the collision of 

Send offprint requests to: P. Guillard, e-mail: two flows of multiphase dusty gas and proposed a model that 

pierre.guillard@ias.u-psud.fr quantifies the gas cooling, dust destruction, H 2 formation and 
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excitation in the postshock medium. H2 gas can form out of gas 
that is shocked to velocities sufficiently low (V s < 200 km s ) 
for dust to survive. Because H2 molecules form on dust grains 
(e.g. Cazaux & Tielens, 2004), dust is a key element in this 
scenario. 

Xu et al. (1999, 2003) reported detection of diffuse far- 
infrared (hereafter FIR) emission from the intergalatic medium 
(hereafter IGM) with the Infrared Space Observatory (ISO), and 
proposed that the dust emission in the shock region would arise 
from dust grains that efficiently cool the X-ray emitting plasma 
via collisions with hot electrons. The FIR emission would then 
trace the structure of the shock, as suggested by Popescu et al. 
(2000) for the case of shocks driven into dusty gas that is accret- 
ing onto clusters of galaxies. However, the poor spatial resolu- 
tion of these observations makes it difficult to separate the dust 
emission associated with star formation (in the neighborhood 
galaxies, or SQ-A) from that really associated with the shock. 
In addition, Guillard et al. (2009) show that the dust contribu- 
tion to the cooling of the hot (~ 5 x 10 6 K) plasma is expected 
to be low, because of efficient thermal sputtering of the grains. If 
we assume that the age of the galaxy collision is ~ 5 x 10 6 yrs, 
grains smaller than 0.1 pm in radius must have been destroyed. 
However, the plasma could still be dusty if, before the shock, a 
significant fraction of the dust mass was in larger grains and/or if 
dust destruction is balanced by mass exchange between the cold 
and the hot gas phases (Guillard et al., 2009). 

The discovery of bright H2 emission in the ridge set new per- 
spectives on the origin of the dust emission in the SQ ridge. In 
the context of our model for the H2 formation in the SQ shock, 
we expect some dust emission from the molecular gas. We use 
Spitzer observations and a model of the dust emission to test this 
expectation. Spitzer observations show that the bright Polycyclic 
Aromatic Hydrocarbons (henceforth PAHs) and mid-IR contin- 
uum emitting regions are spatially correlated with UV emission 
and associated with star-forming regions mainly associated with 
the individual sources in the group (in particular the spiral arm of 
the intruder galaxy NGC 7318b) and with SQ-A (Cluveret al., 
2010; Natale et al., 2010). These IR-bright regions do not corre- 
late with the radio, X-ray, or H2 line emission, which trace the 
shock structure. In this paper we focus on the fainter dust emis- 
sion from the the SQ ridge itself, outside these star-forming re- 
gions. The spatial distribution of the mid-IR emission allows us 
to isolate the dust emission from the shock. UV, visible and near- 
IR observations determine the spectral energy distribution (here- 
after SED) of the radiation eld used as an input to the dust model. 
We consider two limiting cases for the structure of the molecular 
gas, either diffuse, or fragmented into clouds that are optically 
thick to UV light. An updated version of the Desert et al. (1990) 
model is used to compute the dust emission from molecular gas 
in these two cases, and fit the model results to the observed IR 
SED in the SQ shock. The Galactic dust size distribution is taken 
as a reference. 

This paper is organized as follows. Section 2 presents the 
new IR Spitzer observations of SQ, and the UV, optical and 
near-IR ancillary data we use in this paper. The method used to 
perform photometry within the SQ shock region and the results 
are described in sect. 3. Section 4 presents the Spitzer imaging 
and spectroscopy results pertaining to the dust emission in the 
shock structure, emphasizing the PAH properties in the shock. 
The physical framework and inputs of the dust modeling are out- 
lined in sect. 5, and the results are discussed and compared to the 
mid-IR Spitzer observations in sect. 6. In sect. 7 we discuss the 
dust processing in the shock. Then, we present our conclusions 



in sect. 8 and propose new observations to constrain the physical 
structure of the molecular gas in the shock. 

In this paper we assume the distance to the SQ group to be 
94 Mpc (with a Hubble constant of 70 km s _l Mpc~ 2 ) and a 
systemic velocity for the group as a whole of 6 600 km s . At 
this distance, 10 arcsec = 4.56 kpc. 

2. Observations of the Stephan's Quintet shock 

In the following paragraphs we present the new mid-IR 
(sect. 2. 1) observations of SQ, and the ancillary UV (sect. 2.2.1), 
optical (sect. 2.2.2), and near-IR (sect. 2.2.3) data, respectively. 

2. 1 . Spitzer IR imaging and spectroscopy 

Stephan's Quintet has been imaged with the InfraRed Array 
Camera (IRAC, Fazio et al., 2004) at 3.6, 4.5, 5.8, 8 pm, with the 
IRS blue peak-up imager (PUI) at 16 pm, and with the Multiband 
Imaging Photometer for Spitzer (MIPS, Rieke et al., 2004) at 
24 and 70 pm. The 70 pm image has been reported in Xu et al. 
(2008). The pixel sizes are 1.8", 2.45" and 5" at 16, 24 and 
70 pm, respectively. Except for the 70 pm image, we direct the 
reader to Cluveret al. (2010) for a description of the observa- 
tional details and data reduction. The upper right panel of Fig. 1 
shows the 16 pm data from the IRS blue PUI. The bottom left 
and right panels show the MIPS 24 and 70 pm images. 

The SQ shock region was also mapped with the IRS spectro- 
graph (Cluver et al., 2010). The Short-Low (SL) and Long-Low 
(LL) modules of the spectrograph have been used, covering the 
wavelength ranges 5.3 - 14.0 and 14 - 38 pm, with spectral res- 
olution of A/ 5 A — 60 - 127 and 57 - 126, respectively. 

2.2. Ancillary data 

We present the UV, optical and near-IR ancillary data we use to 
determine the radiation eld heating the dust. (sect. 5.2). 

2.2.1. GALEX UV imaging 

The upper left corner of Fig. 1 shows the near-UV (NUV) 
image from the Galaxy Evolution Explorer mission (GALEX, 
Martin et al., 2005). These observations have been reported by 
Xuet al. (2005). The pixel size is 1.5" , the wavelength is A = 
2267 A and the bandwith is AA = 732 A. The far-UV (FUV) 
image (A = 1 516 A, AA = 268 A) is also used but is not shown 
here. 

If one excludes the foreground galaxy NGC 7320 (Sd), most 
of the UV emission is associated with the two spiral mem- 
bers of the group, NGC 7319 and NGC 7318b, and the intra- 
group medium starburst SQ-A. The galaxy-wide shock structure, 
which shows up in Ha, radio and X-ray observations, is barely 
visible on the UV images. Based on the comparison between 
ISO and GALEX data, Xu et al. (2005) concluded that most of 
the UV emission in the ridge is not associated with the large- 
scale shock itself, but with H n regions along the spiral arm of 
the intruder NGC 7318b. 

We use the flux calibration described in the GALEX ob- 
server's guide 1 and in Morrissey et al. (2007). The unit data 
number (DN, or 1 count per second, cps) is equivalent to 108 
and 36yuJy for FUV and NUV, or, equivalently, 1.4 x 10~ 15 and 
2.06 x 10~ 16 erg s cirT 2 A -1 , respectively. We do not apply any 

1 http://galexgi.gsfc.nasa.gov/docs/galex/ 
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Fig. 1. UV and mid-IR observations of Stephan's Quintet. Top-left is a near-UV (2 267 A) GALEX image from Xu et al. (2005), 
top-right is the Spitzer IRS 16 pm peak-up image, bottom-left and bottom-right are the 24 /mi and 70 pm Spitzer MIPS images of 
the SQ group, respectively. For all images, the black contours shows the 4 cr (0.3 MJy sr _1 ) 0-0 S(3) H2 line emission detected over 
the SQ ridge, from Spitzer IRS mapping by Cluver et al. (2010). The magenta contours show the 24 pm emission at a 0.25 MJy sr~' 
level. They are used to identify star forming regions that overlap the H2 contours in the ridge. Photometry is performed within the H2 
contours, excluding or not these star-forming regions. The circles indicate 17" beams where aperture photometry is also performed. 
The yellow circle shows the "ON" position in the SQ shock, centered on a = 22h35'59.8", S = +33°58'16.7". The white circle 
show the "OFF" position we use to subtract the background signal, centered on a = 22h35'57.7", 6 = +33°58'23". 



aperture correction for the extended source UV photometry on 
GALEX images since the PSF full width half maxima (FWHM) 
are 4.9" and 4.2" for the NUV and FUV, respectively. 

2.2.2. HST V-band imaging with WFPC2 

We use V-band data taken with the Wide Field Planetary Camera 
2 (WFPC2) onboard the Hubble Space Telescope (HST). For the 
F569W filter, the data consist of two sets of two images (4x800 s 
exposure in total), each set for a given dithering position. The 
data were first processed by the HST pipeline. After alignment 
of the images with the IRAF data reduction software, a median 
combination of the four images has been taken, followed by the 
removal of the remaining hot pixels by a 3 x 3 pixels median 



filter. We obtained a V-band image similar to that presented in 
Gallagher et al. (2001). 

The photometry has been performed using the flux cali- 
bration given in the header of the images (1 DN= 4.149 x 
10 21 W rrT 2 A" 1 ). The central wavelength of the F569W V- 
band filter is A - 5 644.4 A. 



2.2.3. Near-IR WIRC imaging 

We use near-IR (NIR) data from recent deep observations with 
the Wide field IR Camera (WIRC) on the Palomar 200-inch tele- 
scope (V. Charmandaris, private communication). The WIRC 
images were taken in July, 2009, and processed with the Swarp 
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software 2 . They are 5 - 6 magnitudes deeper than the corre- 
sponding 2MASS images. The zero-point magnitudes are 24.50, 
22.73 and 23.05, so we use flux calibrations of 0.252, 0.830, and 
0.402 /iJy DN 1 for J, H, and K s , respectively. The correspond- 
ing central wavelengths used are 1 .235 + 0.006, 1 .662 ± 0.009, 
2.159 ± 0.01 1 /im. The images show that most of the NIR emis- 
sion in the SQ ridge is associated with the spiral arm of the 
intruder, NGC 7318b. NIR emission associated with the SQ-A 
northern starburst is also detected. 



3. Photometry on IR, optical and UV images 

We describe the method used to perform the photometry on 
Spitzer mid-IR images for the dust emission and on near-IR, op- 
tical and UV images to estimate the radiation field at the position 
of the SQ shock. 



3.1. Method 

3.1 .1 . Dust emission 

The SQ field of view is a crowded region (see sect. 4. 1 for a de- 
scription of the spatial distribution of the dust emission). In or- 
der to isolate the dust emission from the shock itself, we sum the 
signal within regions that are not contaminated by IR-emitting, 
star-forming regions. Since the shock is surrounded by bright 
sources, we do not apply any aperture correction. For compar- 
ison, three different area are used to perform photometric mea- 
surements. 

1. The signal is summed over a circular aperture of 17" in di- 
ameter that is centered on the SQ ridge, in the middle of the 
X-ray emitting shock front. This 17" aperture corresponds to 
the FWHM of the MIPS beam at X = 70 fim. It is marked 
with the yellow circle on Fig. 1 ("ON" position). 

2. The signal is integrated over the SQ ridge within the IT con- 
tours (black line on Fig. 1) that define the shock structure. 
Except SQ-A, star-forming regions within the black H2 emis- 
sion contour are included. They are defined by the 24jt/m 
iso-flux (F24 lim > 0.25 MJy sr~') contour (magenta line in 
Fig. 1). 

3. The signal is summed over the SQ ridge within the H2 con- 
tour but excluding the star forming regions. To do this, we 
subtract from the signal the emission arising from the inter- 
section of the areas within the magenta 24 fim contours and 
black H2 contours. 

To remove the contamination from the halo of the group, we 
choose to estimate the background level in a region close to the 
H2-emitting shock. This "OFF" region is marked with the white 
circle on Fig. 1 . For the 70 fim image, we cannot exclude that the 
flux within the "ON" position is contaminated by the brighter 
sources around it. Thus, we consider our 70 fim flux as an upper 
limit. A more detailed analysis of the MIPS 70 and 160//m will 
be reported in Natale et al. (2010). 

3.1.2. Radiation field 

The radiation field is the integral of the flux over all directions. 
In the U V domain, due to scattering of light, we can assume that 
the radiation field is isotropic and estimate its strength from UV 
photometry at the position of the shock. Thus we measure the 



UV fluxes within the same apertures than for the Spitzer images, 
using the same "OFF" position. 

The UV fluxes are corrected for both foreground galac- 
tic and SQ internal extinctions. For the extinction curve of 
the Galactic diffuse interstellar medium (Ry = 3.1 curve in 
Weingartner & Draine, 2001a), the visible extinction A ( ^ w) = 
0.24 scales to the FUV and NUV Galactic extinctions at GALEX 
wavelengths, respectively Ap^*' = 0.62 and A^*' = 0.72. We 
use the Xu et al. (2005) values for the internal SQ extinction, i.e. 
Ap^ =0.76 and A^» =0.88. 

The optical and near-IR images show that the shock is sur- 
rounded by bright sources, in particular NGC 7318b and nearby 
star-forming regions. Since there is little scattering at these 
wavelengths, the photometry restricted to the shock area is likely 
to underestimate the optical and near IR intensity of the radiation 
field. The choice of the aperture is rather arbitrary in that case 
and only provide a rough approximation of the optical and near- 
IR radiation field. Note that the photometry has been performed 
after removal of Galactic stars (using the DAOPHOT package). 

Based on an average optical extinction of Ay = 0.6 
(Guillard et al., 2009) for the center position in the ridge, and a 
Galactic extinction curve, we apply the following extinction cor- 
rections: Aj = 0.17, A H = 0.1 1, and A Ks = 0.07 for J, H, and K s 
bands, respectively. We use the values of zero-point fluxes from 
Cohen et al. (2003), i.e. 1594 ±28, 1024 + 20, 666.7+ 12.6 Jy for 
J, H, and K s bands. The internal extinction correction applied on 
the HST photometry is Af^w = 0.575. 

3.2. Results 

The quantitative results about the photometry performed on IR 
Spitzer images are gathered in Table 1 . The surface brightnesses 
are given for the three area described above. We estimate from 
the IRS spectrum that the 17 fim S(l) and 8 pm S(4) H2 line emis- 
sion represent respectively 68 % and 62 % of the IRS Peak-Up 
Imager 16//m and the IRAC 8/iin in-band flux within the 17" 
ON aperture centered on the shock, and we correct for this con- 
tamination. The last row of table 1 indicates the fluxes we adopt 
throughout this paper. The error bars are estimated by using two 
different background subtractions close to the "OFF" position 
(8-arcsec shifts in the East- West direction). In addition, we in- 
clude calibration uncertainties (of the order of 5%). 

The results of the UV GALEX photometry are gathered in 
Table 2. The surface brightnesses are given for the three aper- 
tures defined in sect. 3.1.1. The UV fluxes measured over the SQ 
ridge aperture that includes the 24jum-bright (> 0.25 MJy sr -1 ) 
regions are 10 - 25 % higher than the fluxes where these regions 
have been removed from the aperture ("SQ ridge (partial)" row). 
The brightnesses are background-subtracted and corrected for 
foreground galactic extinction and SQ internal extinction. From 
the UV surface brightness values we derive the intensity of the 
standard interstellar radiation field in Habing units. 

The results show that the UV flux in the shock, outside star- 
forming regions, corresponds to an interstellar radiation field 
of average intensity Guv = 1.4 ± 0.2 in Habing units 3 (see 
sect. 5.2.1 for details). The value of Guv is used to character- 
ize the intensity of the non-ionizing radiation field in the shock 
(sect. 5.2). The errorbar on the Guv factor takes into account two 
background subtractions obtained by shifting the OFF position 
by 8" in the East- West direction. This uncertainty does not take 
into account the uncertainty on the optical extinction. Although 



http://astromatic.iap.fr/software/swarp 



The flux of the Habing field (Guv = 1) equals 2.3 x 
1 530 A (Habing, 1968; Mathis et al., 1983). 



10 3 erg s 1 cm 2 at A 
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Table 1. Summary of the mid-IR photometric measurements 3 performed on the Spitzer IRAC, IRS PUI 16 ,um, and MIPS 24, and 
70 /im images (Fig. 1). The background-subtracted surface brightnesses are indicated for three different apertures. The last row (17" 
circular aperture ) shows the values we adopt throughout this paper, and in particular for the SED presented in Fig. 2. 



Target region 


3.6yum 


4.5 yum 


5.1 


/um 


8yum c 


16//m f 


24yum 


70 yum 


SQ ridge b 


23.7 ±0.9 


11.9 ±0.2 


7.3 


±0.2 


4.2 ±0.1 


1.25 ±0.05 


2.25 ±0.12 


< 6.43 


SQ ridge (partial) c 


19.4 ± 1.7 


9.9 ± 1.8 


6.0 


±0.1 


3.35 ±0.11 


1.00 ±0.06 


0.89 ±0.15 


< 5.04 


17" beam d 


17.6 ±0.9 


9.15 ±0.12 


5.29 


±0.15 


3.07 ±0.11 


0.88 ± 0.03 


0.79 ±0.15 


< 3.91 



" Fluxes are in units ofxlO 8 Wm 2 sr 1 . 

b Mid-IR signal integrated over the shock structure, denned by the 4<x S(3) IF rotational line contours (black line on Fig. 1). We exclude the 
SQ-A star-forming region. 

c Mid-IR signal integrated over the S(3) H 2 rotational line contours by excluding overlaps with the (T^^m > 0.25 MJy sr') 24yum magenta 
contours (see Fig. 1). This allows us to partially remove the contribution from star-forming regions to the dust emission. 
d Mid-IR signal integrated over the 17" beam in the center of the SQ shock structure (yellow circle on Fig. 1) 
e The contribution of the 8 //m S(4) Hi line emission to the 8 /jm in-band flux is subtracted. 
' The contribution of the 17/jm S(l) IF line emission to the 16 yum in-band flux is subtracted. 



Table 2. Summary of the UV photometric measurements performed on GALEX FUV and NUV images (see Fig. 1). The flux 
columns indicate surface brightnesses that are corrected for the sky background. 





FUV (A = 


1516 A) 


NUV (A = 


2 267 A) 


Target position 


Flux 


Radiation Field G 


Flux 


Radiation Field G 




xl0~ 7 [ Wnr 2 sr 1 ] 


[Habing units] 


xl0~ 7 [Wnr 2 sr 1 ] 


[Habing units] 


SQ ridge a 


2.82 ± 0.34 


1.54 ± 0.19 


2.74 ±0.39 


1.98 ±0.29 


SQ ridge (partial) b 


2.29 ± 0.35 


1.25 ±0.19 


2.25 ± 0.37 


1.63 ±0.29 


17" ON beam c 


2.49 ± 0.35 


1.36 ±0.19 


2.42 ± 0.39 


1.75 ±0.29 



" UV signal integrated over the shock structure, denned by the 4<x S(3) Hi rotational line contours (black line on Fig. 1). We exclude the SQ-A 
star-forming region. 

* UV signal integrated over the S(3) H2 rotational line contours by excluding overlaps with the (Tu^m > 0.25 MJy sr 1 ) 24 yum magenta contours 
(see Fig. 1). This allows us to partially remove the contribution from star-forming regions to the dust emission. 
UV signal integrated over the 17" beam in the center of the SQ ridge (yellow circle on Fig. 1) 



Table 3. Optical HST V-band and near-infrared WIRC surface brightnesses 3 in the J, H and K s bands for the SQ ridge. 



Band 


V 






J 


H 


K s 


J:H:K S 


Flux 
Dered 


7.2 ± 
12.2 ± 


1.6 

2.9 


8.45 
9.93 


±0.64 
±0.77 


8.54 ± 0.65 
9.44 ± 0.77 


8.07 ± 0.61 
8.60 ± 0.75 


1.048: 1.058 : 1 
1.154 : 1.096 : 1 



" The fluxes are given in units of 10 7 W m 2 sr "Dered" indicates the extinction-corrected brightness, and the "J:H:K S " column shows flux 
ratios between the bands. 



optical spectroscopy shows that there are spatial variations of the 
Ay value in the shock (Ay = 0.1 - 2.5), we do not expect that 
this has an important impact on our estimate of Guv, because it 
is calculated over a large aperture. If the most extincted regions 
have a reasonably small surface filling factor, they will not affect 
Guv significantly. 

The extinction-corrected surface brightnesses measured on 
HST V-band and Palomar WIRC images are gathered in Table 3. 



4. Dust emission from Stephan's Quintet: 
observational results 

This section reports the observational results about the dust 
emission detected in the SQ ridge. We compare the mid-IR SED 



and the relative intensities of the PAH bands in the SQ shock to 
that of the Galactic diffuse interstellar medium (ISM). 



4.1. Spatial distribution of the dust emission from the SQ 
group 

The images in Fig. 1 show that the bright mid-IR emitting re- 
gions are associated with the NGC 7319 galaxy and with UV- 
luminous, star-forming regions (Cluver et al., 2010). The mid-IR 
emission does not correlate with the FF, X-ray and radio emis- 
sions. Interestingly, these star-forming regions are outside the 
galactic disks of NGC 7318a and b, suggesting that a significant 
amount of gas has been displaced from these galaxies by tidal in- 
teractions. If we exclude the foreground galaxy NGC 7320, the 
bright mid-IR emission comes from the NGC 7319 galaxy, the 
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SEDs of the SQ shock and Galactic diffuse ISM 



IRAC Bands 
I I I I 



i-B-i 



i-B-i 
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Stephan Quintet shock A 
Galactic diffuse ISM (x 1/200) □ 



10 100 
Rest-frame wavelength [/xm] 

Fig. 2. The SED of the dust emission from the SQ shock seen 
by Spitzer (table 1), compared with that of the Galactic dif- 
fuse emission observed with Spitzer IRAC, IRAS and DIRBE 
towards the line-of-sight centered on the Galactic coordinates 
(28.6, +0.8). The contribution of the gas lines is subtracted. The 
Galactic SED is scaled down by a factor 200. Between 3-8 pm, 
the SQ flux is much brighter because of the stellar contribution. 
Horizontal bars indicate the filters bandwiths. 



SQ-A starburst region, and an arc-like feature to the east of the 
intruding galaxy NGC 7318b. This arc structure is clearly seen 
on the UV, IRAC, and 16ji/m images, and could be associated 
with NGC 7318b's spiral arm. The 24 and 70yum images show 
a bright and extended emission at the southern tip of the ridge 
which may be associated with enhanced star formation in this 
arc feature. 

Dust is detected in the SQ shock region, at a distance of 
10-20 kpc of the nearest surrounding galaxies. In the following, 
we focus on this faint emission observed within the SQ ridge, 
outside star-forming regions. Note that on the 16 pm IRS Blue 
Peak-Up image, the emission in the ridge correlates well with 
the H2 emission, which is due to the contamination of the H2 
17 pm S(l) line to the in-band flux (sect. 3.2). 

Using a combination of the Ha and 24 pm luminosities, or 
the 1.1pm PAH emission, Cluveretal. (2010) find an upper 
limit of < 0.08 M yr" 1 on the star formation rate in the shock, 
as compared with ~ 1.25 M Q yr" 1 in SQ-A. This suggests that 
the star formation is being depressed in the shock region. 

4.2. Dust emission from the SQ shock, outside star-forming 
regions 

The images and photometric measurements show that thermal 
dust emission is detected in the SQ shock structure, outside star- 
forming regions. The SED of the IR emission from the center of 
the shock is shown on Fig. 2. Among the SEDs listed in table 1, 
the figure displays that of the last row (17" beam, see sect. 3). 

For comparison with Galactic data, the figure includes a 
12 pm flux computed by integrating the IRS spectrum over 
the IRAS 12//m filter bandpass after subtraction of the gas 
lines. The 12/24 pm flux density ratio in the SQ shock is 
T^ Q \\2pm)llf Q \2Apm) = 0.97 ± 0.15, which is re- 
markably close to the value for the Galactic diffuse ISM: 

0.07. The Galactic 



SED is that measured towards the line-of-sight centered 
on the Galactic coordinates (28.6, +0.8), observed with the 
ISOCAM-CVF 4 and Spitzer IRAC (Flagey et al., 2006). Fig. 2 
also shows the IRAS 12, 25, 60 and 100 pm fluxes for 
the Galactic diffuse ISM emission measured on IRIS images 
(Miville-Deschenes & Lagache, 2005). The IRAS SED is ex- 
tended to 140 and 240 pm using DIRBE (Hauser et al., 1998) 
color ratios at the position of the line of sight. The total col- 
umn density for this line of sight is estimated to be A^h = 
2 x 10 22 cm" 2 , and the mean radiation field a few Habing units 
(Guv ~ 3). In the figure, the Galactic SED is scaled down by a 
factor 200. 

We propose that the faint dust emission we have isolated in 
the center of the SQ shock is diffuse emission associated with 
the shocked molecular gas present in the ridge. An second pos- 
sibility is that the dust is associated with the Hi gas and a third 
one that it is produced by unresolved star-forming regions. We 
favor the first interpretation for three reasons. 

- Looking at the Hi data of Williams et al. (2002), the two 
outer contours in their figure 9, 0.6 and 1 x 10 20 cm" 2 , inter- 
sect our 17" aperture used for dust photometry. Since the an- 
gular resolution of these Hi observations is 20", the Hi emis- 
sion is likely to be contaminated by beam smearing of the 
brighter emission to the north of our aperture. On the south- 
ern side of our aperture, in the ridge, Hi is undetected and 
Nui < 5.8 x 10 19 cm" 2 . This upper limit is smaller than the 
column density of warm H 2 (2 x 10 20 cm" 2 ) derived from 
the Spitzer H2 fluxes. Since the warm H2 column density 
is a lower limit on the total H2 column density, Hi gas ac- 
counts for a minor fraction of the gas column density in the 
ridge. Therefore the dust emission cannot be mostly associ- 
ated with the Hi gas in the shock region. 

- HST observations show that there are very few star clus- 
ters in the center of the ridge. Most of them are associ- 
ated with NGC 7319, the tidal debris of NGC 7318a/b, 
and the SQ-A intragroup starburst region (Gallagher et al., 
2001). We find three candidates for star clusters (My = 
-12.14, -10.01, -9.65) within our 17" (~ 7.8 kpc) aper- 
ture in the center of the shock. These three clusters produce 
a V-band flux of 5.2 x 10" 9 W m" 2 sr" 1 , which is 2 orders of 
magnitude lower than the V-band surface brightness that we 
derived from HS T observations in the center of the ridge. 

- The 12/24yum flux ratio is remarkably similar to that of the 
Galactic diffuse ISM. The average column density of warm 
molecular gas in the SQ shock derived from H2 observations 
is Nh = 2 x 10 20 cm" 2 , a factor 100 smaller than that of the 
Galactic line of sight, in agreement with the scaling factor 
used to match the fluxes in Fig. 2. 



4.3. Mid-IR spectroscopy: characterization of the dust 
emission from the SQ shock 

4.3.1 . Spitzer IRS spectrum in the center of the shock 

The top panel of Fig. 3 presents a low-resolution Spitzer IRS 
spectrum extracted from a central position of the SQ ridge. This 
spectrum was obtained by Cluver et al. (2010), averaging all the 
spectra observed within a 274.2 arcsec 2 (~ 18" x 15") rectan- 
gular aperture around the center of the ON position. These new 
data has much higher sensitivity and better flux calibration than 
the first observations reported by Appleton et al. (2006). This 



r! MW) (12pm)/T^ W) (25pm) = 0.95 



4 Camera equiped with a Circular Variable Filter (CVF) on-board the 

Infrared Space Observatory (ISO) 
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Spitzer IRS spectrum extrocted in the SQ ridge 
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Fig. 3. Spitzer IRS mid-IR spectra extracted over a ~ 18" x 15" area in the center of the SQ shock structure (RA 22:35:59.9, DEC 
+33:58:16.8). This aperture corresponds approximately to the "ON" beam shown in Fig. 1. Top: full spectrum. Bottom: gas lines 
are removed to highlight PAH dust features and thermal dust continuum. Photometric measurements (red triangles) performed on 
mid-IR images (same as in Table 1 and Fig. 2) are overplotted. Horizontal bars indicate the filters bandwiths. 



PAH spectrum: SQ shock and diffuse Galactic medium 
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Fig. 4. Spitzer IRS PAH emission spectrum extracted in the center of the SQ ridge (black solid line), compared with the ISOCAM- 
CVF spectrum of the diffuse galactic medium (blue dashed line, with flux density labeled on the right). The IRS spectrum is 
smoothed to a resolution H — A/6 A — 24 — > 51, comparable with the ISO- CVF resolution CR = 35 — > 45). Gas lines are removed in 
both spectra. Note that the flux scales are different for each spectrum (labels on the left for SQ, and on the right for the Galaxy). 
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spectrum confirms the first results discussed in Appleton et al. 
(2006). It shows bright S(0) to S(5) H 2 rotational lines and for- 
bidden atomic lines, with a remarkable 34.8 fim [Sin] feature. 
Here we focus on the dust emission. 

The IRS spectrum shows that PAH and thermal dust emission 
is detected from the center of the SQ ridge, outside star-forming 
regions. A weak thermal continuum is also visible from 20 to 
38 fim. The ratio between the fluxes of the 17 fim S(l) H 2 line 
and the 7.7 fim PAH feature is 1.2 + 0.2 in the shock, which is 
about two orders of magnitude higher than the value observed in 
star-forming galaxies (Roussel et al., 2007). 

4.3.2. PAH emission from the SQ shock 

The 7.7, 11.3 and 17 /mi Aromatic Infrared Bands (AIBs), at- 
tributed to PAHs, are detected (although weak) in the new 
Cluveret al. (2010) observations at the center of the SQ shock. 
We use the PAHFit IDL tool 5 (Smith et al., 2007) for decompos- 
ing the full 5-38 fim IRS spectrum into contributions from PAH 
features, thermal dust continuum, starlight and gas lines. We do 
not include any extinction in the fit. The result of the fit on the 
full 5-38 fim spectrum is shown in Appendix A, Fig. A. 1 . This 
spectral decomposition allows us to remove the gas lines and 
extract a "pure" dust spectrum, shown on the bottom panel of 
Fig. 3. PAHFit is run one more time on this spectrum, which al- 
lows us to measure accurately the fluxes of the AIBs. The line 
strengths of the PAH emission features and their ratios are gath- 
ered in Tables 4 and 5, respectively. 

In Fig. 4 we compare the PAH spectrum from the SQ shock 
with the ISOCAM-CVF spectrum of the diffuse Galactic emis- 
sion from Flagey et al. (2006), between 5 and 20 fim. The main 
differences between the two spectra are the following. First, 
the ratio between the flux in the band at 7.7 /mi and the band 
at 11.3 /mi (henceforth ^7.7/11.3) is a factor ~ 2 lower for SQ 

(R^y n3 = 1.36 ± 0.07, see table 5) than for the diffuse Galactic 

light (4j)?i.3 = 2 - 93 ± °- 08 )- The #7.7/11.3 value for the SQ 
shock is comparable to that measured for AGN of the SINGS 6 
sample (Smith et al., 2007). Second, the 6.2 fim AIB is absent 
in the SQ spectrum. Note that the rise of the SQ spectrum at 
A < 6 fim is due to the stellar component. The 17 fim complex is 
prominent but the 16.4pm feature is not seen on top of it. This 
may not be significant because, to our knowledge, this feature is 
only prominent in the NGC 7023 Galactic PDR 7 (Sellgren et al., 
2007). This feature is also observed in the star-forming galaxy 
NGC 7331 (Smith et al., 2004), but with a higher signal to noise 
spectrum than the one we have for SQ. 

The enhancement of the 1 1.3 fim AIB compared to the 6.2, 
7.7 and 8.6 fim features has been observed on galactic scales in 
elliptical galaxies (e.g. Kaneda et al., 2005), and in active galac- 
tic nuclei (e.g. Smith et al., 2007). This enhancement has also 
been discussed in the Galactic diffuse medium (e.g. Flagey et al., 
2006) or at small scales, in PDR interfaces (e.g. Rapacioli et al., 
2005; Compiegne et al., 2007). 

The PAH emission spectrum depends on the size distribu- 
tion, hydrogenation and ionization states. The Ru/u.3 PAH band 
ratio depends mainly on the PAH ionization state. Theoretical 
(Langhoff, 1996; Bakes et al., 2001a,b; Draine & Li, 2001; 
Bauschlicher, 2002) and experimental (e.g. Szczepanski & Vala, 
1993) studies show that neutral PAHs have lower ^7.7/11.3^ 

5 available on http://tir.astro.utoledo.edu/jdsmith/pahfit.php 

6 Spitzer Infrared Nearby Galaxies Survey, http://sings.stsci.edu/. 
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than charged ones (both anions and cations). The charge state 
of PAHs is mainly determined by the ionization parameter 
Guv Vr / "e where Guv is the integrated far ultraviolet (6 - 
13.6 eV) radiation field expressed in units of the Habing radi- 
ation field, T is the electron temperature and n e is the electronic 
density. This parameter translates the balance between photoion- 
isation and recombination rates of electrons (Bakes & Tielens, 
1994; Weingartner & Draine, 2001c). Flagey et al. (2006) quan- 
tified ^7.7/1 up,,, as a function of the ionization parameter 

Guv /«e and the PAH mean size. We use their calculations 
to discuss the ionization state of the PAHs in the SQ shock. Our 
PAHFit decomposition yields Rj.j/uium = 1.36 + 0.07, which 
translates into Guv ^lT ln e < 30 K 1/2 cm 3 . Assuming a warm 
molecular gas temperature of ~ 150 K and Guv = 1, we find 
that « e > 0.4 cm -3 . This lower limit is one order of magni- 
tude higher than the electronic densities inferred from observa- 
tions and modeling of the ionization of cold neutral gas in the 
Solar neighbourhood (Weingartner & Draine, 2001b). If this in- 
terpretation and diagnostic applies, a higher ionizing flux from 
cosmic -rays or X-rays would be required to maintain such a high 
electron density. 

The /?6.2/7.7 PAH flux ratio depends on the size of the emit- 
ting PAHs (Draine & Li, 2001). The non-detection of the 6.2 fim 
band sets a low upper limit of R^.ip.i < 0.2 on the PAH 
strength ratio, that suggests preferentially large PAHs. Although 
the method used by Draine & Li (2001) to measure the PAH line 
strengths is different from ours, we find that both ^7.7/11.3 and 
#6.2/7.7 flux ratios in SQ can be explained by large (with a num- 
ber of carbon atoms N c > 300) and neutral PAHs in CNM con- 
ditions, excited with a Guv — 1 interstellar radiation field (see 
Fig. 17 of Draine & Li, 2001). 



5. Modeling dust emission 

We present the physical framework and inputs of our model- 
ing of the emission from dust associated with the molecular gas. 
Section 5.1 and 5.2 present the codes and the radiation field we 
use for our calculation of the dust emission. 

5.1. The DUSTEM (Dust Emission) code 

We use an updated version of the Desert et al. (1990) model, 
the DUSTEM code, to compute the dust emission. The modi- 
fications 8 implemented to the Desert et al. model are described 
in Compiegne et al. (2008). The dust properties, the dust-to-gas 
mass ratio, and the incident radiation field being given, the code 
calculates the dust SED vS v in units of erg s _1 H" 1 , for each dust 
grain species, as a function of the wavelength. 

We use the diffuse Galactic ISM size distribution (a power- 
law n(a) oc gt 3 - 5 , Mathis et al., 1977), and dust-to-gas mass ra- 
tio infered from the fitting of the SED and extinction curve of 
the diffuse ISM (Compiegne et al., 2008). The DUSTEM model 
includes a mixture of three populations of dust grains with in- 
creasing sizes: 

- Polycyclic Aromatic Hydrocarbons (PAHs) of radius a - 
0.4 - 1.2 nm, responsible for the Aromatic Infrared Bands 
(AIBs) and the FUV non-linear rise in the extinction curve. 



8 In particular, the absorption cross sections of the PAHs (with addi- 
tion of new Aromatic Infrared Bands, AIBs) and the Big Grains (BGs), 
as well as the heat capacities (graphite, PAH C-H, silicate and amor- 
phous carbon) have been updated. 
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Table 4. Fluxes of the PAH bands measured with the PAHFit IDL tool on the Spitzer IRS spectrum of the center of the SQ ridge and 
the ISO-CVF spectrum of the Galactic light. 



Integrated fluxes (see notes for units) 


PAH band 


6.2 pm 


7.7^m c 11.3//m c 12.6yum c 


17 /xm c 


SQ ridge (IRS) 11 


< 0.75 


3.8 ±0.1 2.8 ±0.1 1.2 ±0.2 


0.93 ± 0.04 


Galactic (CVF) b 


5.19 ±0.09 


19.05 ± 0.03 6.49 ± 0.17 2.90 ± 0.05 


No data 



" Fluxes are in units of x 10~ 9 [W irT 2 sr~'] 
* Fluxes are in units of xl0~ 7 [W nT 2 sr~'] 

c In the case of PAH blended complexes, we sum the fluxes of the different lorentzian components that contribute to the feature. 



Table 5. Flux ratios between the PAH bands for the Spitzer IRS spectrum of the center of the SQ ridge and the ISO-CVF spectrum 
of the diffuse Galactic light. 



Flux ratios 




Rb.ip.l 


#7.7/11.3 ^ 1 1.3/ 12.6 


^11.3/17 


SQ ridge (IRS) 


< 0.2 


1.36 ±0.07 2.3 ±0.2 


3.01 ±0.16 


Galactic (CVF) 


0.272 ± 0.005 


2.93 ± 0.08 2.24 ± 0.07 


No data 



- Very Small Grains (VSGs, a — 1-4 nm), which are carbona- 
ceous (graphitic) nanoparticles producing the mid-IR contin- 
uum emission and the extinction bump at 2175 A. 

- Big Grains (BGs, a = 4-110 nm) of silicates with car- 
bonaceous mantles or inclusions, which account for the far 
IR emission and the l/A rise at visible and near-IR wave- 
lengths. 

5.2. Radiation field 

In this section we model the radiation field used to compute the 
dust emission. Stellar radiation is coming from the surround- 
ing galaxies and/or stars in the ridge (Gallagher et al., 2001). 
The presence of ionizing radiation in the SQ ridge is also in- 
dicated by the lack of Hi gas in the H2-bright shock structure 
(Sulentic et al., 2001), and by emission from ionized gas lines. 
Optical line emission diagnostics suggest that shocks are respon- 
sible for hydrogen ionization in the ridge (Xu et al., 2003). 

Therefore, we consider that the SED of the input radi- 
ation field consists of two components: a stellar component 
(sect. 5.2.1), and photo-ionizing radiation from gas shocked at 
high-velocities (sect. 5.2.2). 

5.2.1 . The interstellar radiation field (ISRF) 

The photometry performed on NUV and FUV GALEX images 
shows that UV flux in the shock, outside star-forming regions, 
corresponds to an interstellar radiation field of intensity Guv = 
1 .4 + 0.2 (see Table 2). What is the origin to this U V radiation in 
the shock? We estimate the UV emission in the ridge that comes 
from Hn regions associated with star formation in the surround- 
ing sources, i.e. NGC 7318a/b, the star-forming region SQ-A 
at the northern tip of the shock structure, and NGC 7319. The 
extinction-corrected FUV luminosities of these main sources 
that surround the ridge are L FUV (NGC7318ab) = 8.1 x 1O 9 L , 
^fuv(SQ-A) = 6.9 x 1O 9 L , and L FUV (NGC7319) = 6.9 x 
1O 9 L (Xuetal., 2005). Assuming that their distances to the 
center of the ridge are 10, 20 and 25 kpc, respectively, we find 
that the total UV flux in the ridge is 3.5 x 10~ 3 erg s _1 cm -2 , 



i.e. Guv - 1-5 in Habing units. Most of the UV flux comes 
from NGC 7318a/b. Although this simple calculation does not 
take into account the extinction between the UV sources and the 
ridge, it shows that this estimate of the UV field is roughly con- 
sistent with our estimate from UV observations. The very small 
number (of the order of unity) of star clusters that lie within the 
center of the SQ shock (Gallagher et al., 2001) shows that their 
contribution to the UV field is small. 

The SED of the stellar component of the radiation field is 
shown on Fig. 5 (green dashed line). We use the Mathis et al. 
(1983) ISRF, scaled by a factor 1.4 to fit the FUV and NUV 
GALEX photometry. The HST V-band and WIRC near-IR pho- 
tometric measurements (see sect. 2.2.3 for observational details) 
are overlaid on Fig. 5 (see Table 3). The colors of the J, H and K s 
fluxes do not exactly match that of the (Mathis et al., 1983) radi- 
ation field, but we remind that these measurements are uncertain 
(because at these wavelengths the radiation field is anisotropic, 
see sect. 3.1.2). This effect does not affect significantly our dust 
modeling since the near-IR part of the radiation field has a small 
impact on the dust emission. 

5.2.2. Photo-ionizing radiation field: shock and precursor 
components 

In the shock sub-region, the [Oi]6300A and [Nn]6584A to 
Ha line ratios are ~ 0.7 and ~ 0.3, respectively, which is 
evidence of shock ionization (Xu et al., 2003). Therefore we 
model the SED of the ionizing field using emission from 
a radiative shock. In the following we constrain the shock 
parameters we use to model the emission from the ion- 
ized gas. The mid-IR IRS spectrum extracted in the core 
of the shock (Fig. 3) shows fine-structure line emission 
from [Nen]/l 12.81 /an, [Nem]i 15.56//m, [Fen]^ 25.99 pm, 
[S m]A 33.48 pm, and [Si n]A 34.82 pm. The spatial distribution 
of this emission and mid-IR line diagnostics of the excitation 
mechanisms are discussed in Cluveret al. (2010). We summa- 
rize here the main results that are relevant to constrain the shock 
parameters. 

The comparison between the high value of 
^[Neiijia.s/^Neiirjis.e ~ 3.55 and shock models (Hartigan et al., 
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Fig. 5. The 2-component SED of the radiation field (stellar + shock) used as input for dust models to calculate the dust emission 
from molecular gas. The radiation field is the sum of an ISRF of intensity Guv = 1.4, scaled to match the observed UV flux in the 
shock, and the photo-ionizing spectrum of a V s = 100 km s~' shock, scaled to match the Ha emission in the SQ ridge. The preshock 
gas density is «h = 1 cirT 3 and the preshock magnetic field strength is Bo = 1 /^G. 



1987; Allen et al., 2008) allows to firmly constrain the range of 
shock velocities to ~ 80 - 200 km s -1 (Cluveret al., 2010). In 
the case of clumpy gas, if we consider the line emission from 
the shock only (discarding line emission from the pre-shock gas 
ionized by emission from the shocked gas), the upper limit on 
shock velocities is a little higher (300 instead of 200 km s~'). 
In addition, comparison of the [S ra]/l 33.5 jt/m/[S m]A 18.7 fim 
line ratios with these shock models constrains the pre-shock 
gas density to be «h <; 10 cm -3 . Therefore, we adopt in the 
following a shock velocity of V s = 100 km s -1 and a pre-shock 
density of «h = 1 cm -3 . The photo-ionizing emission spectrum 
from shocked gas is taken from the library of the Mappings 
III shock code 9 (Allen et al., 2008). We normalize the shock 
spectrum to the observed Ha flux 10 . After this normalization, 
the SED of the ionizing shock emission is weakly sensitive to 
the two main shock parameters, the shock velocity and the gas 
density. 

Fig. 5 shows the SED of the radiation field we use as input 
for the DUSTEM code to compute the outcoming dust emission 
from the molecular and ionized gas. The black solid line is the 
sum of the two contributions (stellar + shock) to the radiation 
field. The shock spectrum itself is mostly composed of ther- 
mal bremsstrahlung (free-free) continuum and resonance lines 
arising from many different elements and ionic stages. It also 
shows a prominent low-temperature bound-free continuum of 
hydrogen, produced in the cool, partially-ionized zone of the re- 
combination region of the shock, and the strong hydrogen two- 
photon continuum produced mostly by the down-conversion of 
Lya photons trapped in this same region of the shock structure. 

9 http://cdsweb.u-strasbg.fr/~allen/mappings.page1 .html 

10 r Ha = 4.1 x 10 -9 W nr 2 sr 1 fromXu et al. (2003). 



Also present, though to a much weaker scale, is the bound-free 
continuum arising from the heavier elements, with the helium 
continuum the most obvious. 

In the case of fast shocks (V s > 100 km s ), the shocked 
medium emits UV radiation that ionizes the pre-shock medium 
before it is shocked. A so-called radiative precursor propagates 
ahead of the shock, with an ionization front velocity that exceeds 
that of the shock. The contribution of the precursor itself is in- 
dicated on Fig. 5 by the red dotted line. The green dashed line 
shows the contribution of the Guv = 1.4 ISRF. The UV GALEX, 
Ha line, and near-IR WIRC fluxes in the center of the shock are 
indicated. 

The contribution of the precursor to the total Hn column 
density and radiative flux depends on the dumpiness of the pre- 
shock medium. This contribution is negligible if the molecu- 
lar gas is clumpy, i.e. fragmented into dense («h > 10 3 cm -3 ) 
clouds that have a small volume filling factor, because most of 
the ionizing photons do not interact with neutral gas but with the 
volume-filling, hot plasma that is optically thin to ionizing radia- 
tion. However, diffuse H2 gas is expected to have a much higher 
volume filling factor. Our modeling of the dust emission from 
the ionized gas takes into account the precursor contribution, 
weighted by the volume filling factor of the clumpy molecular 
gas. 

6. On the structure of the molecular gas 

The dust emission depends on the optical thickness of the molec- 
ular gas to UV radiation. Is the molecular gas diffuse or frag- 
mented in optically thick clumps? In section 6. 1 we describe the 
assumptions we make regarding the structure of the molecular 
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Fig. 6. Sketch of the two limiting physical states of the multi- 
phase molecular gas we are considering in the framework of our 
modeling of dust emission. (Left) Dust is associated with diffuse 
molecular gas that is broken into fragments, filaments or sheets, 
penetrated by UV radiation. (Right) The dusty molecular gas is 
in clouds that are optically thick to UV light. 

gas. Then we present the results of our modeling of the dust 
emission in the SQ shock. We discuss the dust emission for the 
two physical structures of the H2 gas presented in sect. 6.1: dif- 
fuse (sect. 6.2) and clumpy molecular gas (sect. 6.3). Both mod- 
els include the contribution of the ionized gas (H 11) to the dust 
emission (this calculation is detailed in sect. 6.4). This detailed 
modeling is tted to mid-IR observations, and used to investigate 
the inuence of the structure of the molecular gas on the FIR dust 
SED. 

6. 1 . Two limiting cases for the structure of the molecular gas 

The spectral energy distribution of the dust emission depends 
also on the structure of the molecular gas, and in particular on the 
optical depth of the clouds. The physical structure of the molec- 
ular gas in the SQ shock is an open question. In the following 
we explore these two cases. 

1. The molecular gas is diffuse as the gas observed in 
the solar neighborhood through UV spectroscopy (e.g. 
Rachford et al., 2002) and across the Galaxy through obser- 
vations of mid-IR H2 line emission (Falgarone et al., 2005). 
In this case, we assume that the molecular gas is optically 
thin to UV radiation. 

2. The molecular gas is within clouds that are optically thick to 
UV photons, as those observed in star-forming regions. 

Fig. 6 illustrates, in a simplistic manner, the two physical states 
(diffuse or clumpy) of the molecular gas we are considering. In 
both cases, the molecular gas is embedded within Hn gas and 
shock-heated, X-ray emiting plasma. 

In order to compute the dust emission from optically thick 
(to UV photons) clouds (Fig. 6, right panel), we use the Meudon 
PDR (Photon Dominated Region) code (Le Petit et al., 2006) 
to compute the radiative transfer through the cloud. This 1- 
dimensional, steady-state model considers a stationary plane- 
parallel slab of gas and dust, illuminated by UV radiation. The 
radiation field output I v of the Meudon PDR code is used as input 



to the DUSTEM program to compute the spectral energy distri- 
bution (SED) of the dust as a function of the optical depth into 
the cloud. This 2-step process is iterated to take into account 
dust heating by the dust IR emission. Usually 4-5 iterations are 
needed so that the radiation field converges. 

6.2. Dust emission from diffuse molecular gas 

The DUSTEM code is used to model the dust emission from dif- 
fuse molecular gas (sect. 5.1), penetrated by UV photons. The 
cloud is illuminated by a composite field, made up of a stellar 
component (ISRF of intensity Guv = 1 .4 Habing units) and the 
ionizing radiation from the shocked gas (see sect. 5.2 and Fig. 5). 
The dust properties are those of sect. 5.1. The abundances and 
size distribution of the dust grains are Galactic. The dust emis- 
sion from the ionized gas is included in the model SED. It is 
calculated separately with the DUSTEM code and added to the 
contribution of the cloud (see sect. 6.4). 

The top panel of Fig. 7 shows the 0.1 - 300 /im SED result- 
ing from the model (black solid line), so that both the UV ra- 
diation field and IR dust emission can be seen. The dashed line 
(in the bottom right corner of the plot) shows the dust emission 
from the Hn gas. The three other broken curves show the con- 
tributions of the different populations of dust grains, i.e. PAHs, 
VSGs and BGs, respectively. The GALEX, WIRC and Spitzer 
fluxes are indicated for comparison (see sect. 3, Tables 1, 2 and 
3). The 7.7, 11.3 and 17 pm points are the peak values of the PAH 
bands detected in the IRS spectrum (sect. 4.3). Other points come 
from imaging broadband photometry (GALEX, HST, IRAC and 
MIPS) measurements performed in the center of the SQ ridge 
(over an aperture of 17" in diameter, see sect. 3). 

The 11.3, 16 and 24 pm data points are used to fit the model 
SED. As it may be a lower limit, the 70pm point has not been 
included in the fit. The gas column density is the only free pa- 
rameter to fit the data. The model SED is obtained by multiply- 
ing the emissivity output of the DUSTEM code by the column 
density of warm molecular gas that is determined from fitting 
the SED to the mid-IR data. We find that the column density that 
best fits the data is N H = 1.8 ± 0.5 x 10 20 cnT 2 for a Galactic 
dust-to-gas mass ratio. This column density obtained from the 
model is close to the column density derived from the H2 line 
fluxes (Nu - 2.2 x 10 20 cnT 2 ) within this aperture. This column 
density is obtained by fitting the rotational H2 line fluxes with C- 
shocks models as described in Guillard et al. (2009). This sug- 
gests that the column density of cold (T < 50 K) molecular gas 
is much smaller than that of the warm H2. However, this may not 
be the case because the dust-to-gas mass ratio is possibly lower 
than the Galactic value we have assumed here, and the molecu- 
lar gas may be in clumps optically thick to U V photons (see next 
section). 

6.3. Dust emission from clumpy molecular gas 

We model the dust emission associated with molecular gas frag- 
mented into clumps that are optically thick to UV radiation. Here 
the DUSTEM and Meudon PDR codes (see sect. 5.1) are com- 
bined to calculate the emission from a molecular cloud of to- 
tal Ay = 3. The column density of the cloud is thus JVg = 
6.9 x 10 21 crrT 2 . The choice of the Ay value is arbitrary. We 
check that the model SED in the mid-IR domain does not de- 
pend much on the total Ay because the mid-IR emission comes 
mainly from the surface of the cloud. The cloud is illuminated 
by a composite field (stellar + shock), made up of the ISRF of 
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Fig. 7. SED of the dust emission associated with diffuse (top) or clumpy (bottom) molecular gas. The cloud is exposed to a radiation 
field consisting of a mixture (stellar + shock) of the Mathis et al. (1977) ISRF (scaled by a factor Guv = 1-4) and the ionizing 
emission spectrum of shock-heated gas for a shock velocity of 100 km s _1 . (top) The dashed line shows the contribution of the 
ionized gas (see text for details). The other broken lines show the different contributions of the three populations of grains (PAHs, 
VSGs, and BGs). We do not show them on the bottom figure for clarity. Assuming a Galactic dust-to-gas mass ratio, a column 
density of Nh = 1.8 + 0.5 x 10 20 cirT 2 is needed to fit the Spitzer data, (bottom) The black solid line is the sum of the incoming and 
outcoming emission from a slab of gas of total Ay = 3, calculated with the Meudon PDR and DUSTEM codes. To fit the Spitzer 
data, the model IR emission is scaled by a factor / = 0.03, which represents the effective surface filling factor of the molecular gas. 
The observed IRAC and MIPS fluxes correspond to the entries of Table 1. The green points are the IRS fluxes at the peak of the 7.7, 
11.3 and 17 fxm PAH bands. 
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intensity Guv = 1.4 (Habing units) and the ionizing radiation 
from the shocked gas (see sect. 5.2 and Fig. 5). 

Fig. 7 shows the sum of the outcoming radiation (cloud + 
ionized gas) emission spectrum, plus the incoming radiation, 
from the UV to the FIR (0.1 - 300 /im). The dashed line shows 
the emission from the ionized gas. The UV, and IR photomet- 
ric measurements are overlaid for comparison between data and 
models. 

From UV to FIR wavelengths, the total SED (solid black 
line) consists of the ISRF, including the free-free, free-bound 
and resonance line emission from the shocked gas, and the dust 
emission including the AIBs, the emission from the VSGs and 
the grey body of the BGs at long wavelengths. The spectrum also 
shows some fine-structure IR lines, indicated on the spectrum, 
e.g. 18.7 pm [S m], 34.8 pm [Si n], and 157.7 pm [C n]. To fit the 
data, the model IR emission from the cloud is scaled by a factor 
of / = 0.03, which represents an effective surface filling factor 
for the molecular gas. The low value of / is consistent with the 
low value of the average optical extinction in the shock region. 
For clumps of Ay = 3, the average column density of the cloud 
phase is thus (A^h) -fx N£ - 2.1 x 10 20 cirr 2 . Our specific 
choice of Ay = 3 leads to a value of {Nh) that is consistent 
with the column density of the warm H2 gas seen by Spitzer. For 
larger values of Ay, the total column density of molecular gas 
will be larger than the column density of warm H2 inferred from 
rotational lines. 



6.4. Dust emission from shock- and precursor-ionized gas 

We detail how we calculate the emission from dust associated 
with the ionized gas. This contribution is added to both models 
(diffuse and clumpy). We use the input radiation field of Fig. 5, 
including the ionizing part of the spectrum (A < 912 A), to com- 
pute the dust emissivity with the DUSTEM code. The column 
density of ionized gas has been constrained by the integrated 
flux of the 12.8 /im [Nen] fine-structure line, ^[Neii], measured 
on the IRS spectrum. The ^[Neinj/^Neii] line flux ratio is low 11 , 
which implies that Ne 11 is the main ionization state of Ne in the 
shock. Therefore, the column density of ionized gas can be ex- 
pressed as: 



^(ion.gas) = 2 3 x 1Ql9 



r [N eii] [Wm^sf 1 ] 10 cm 



1.17 x 10" 



n e 



(1) 



We use the relation between the [Nen] line intensity and the 
emission measure given in Ho & Keto (2007) and we assume 
that the electronic density of the 10 4 K gas is n e = 10 cm -3 . 
This value is consistent with the SQ postshock pressure (~ 
2 x 10 5 K cm 3 ). The column density derived from Eq. 1 is 
comparable to the H 11 column density derived from the shock 
model (A^h(Hii)= 1.2 x 10 19 cm -2 for a shock + precursor 
model at a shock velocity of 100 km s and a preshock den- 
sity «h = 1 cirr 3 ). 

The dust emission associated with this amount of ionized 
gas is the black dashed line on the top plot of Fig. 7. The column 
density of the ionized gas being one order of magnitude smaller 
than that of the warm H2, the contribution of the ionized gas to 
the IR dust emission is negligible. 



[Nelll] 



ell] 



0.14 ± 0.04 in the center of the SQ ridge, see 



Table 3 in Cluver et al. (2010). 



6.5. Comparison between models and degeneracies 

Fig. 8 shows the comparison between the result of the two mod- 
els presented above. We focus on the dust emission, so the gas 
lines and the stellar continuum from the ISRF are removed from 
the spectra. The FIR SEDs are different between the two models. 
The FIR peak brightness of the SED of the diffuse molecular gas 
is brighter (by a factor of 1.7) and shifted towards shorter wave- 
lengths compared with the clumpy model. This difference can 
be easily explained. When the molecular gas is clumpy, the dust 
grains see, on average, an attenuated radiation field, and are thus 
colder than in the diffuse model. For a given column density of 
matter, the FIR brightness is also fainter. 

Present observations do not allow us to decide whether the 
gas is diffuse or fragmented into optically thick clumps. The spa- 
tial resolution of Spitzer at long wavelengths {A > 70 pm) is not 
high enough to obtain accurate photometric measurements. FIR 
observations with the Herschel Space Telescope will provide 
the sensitivity and angular resolution needed to test our mod- 
els. However, it may not be so straightforward to conclude on 
the structure of the molecular gas because we assume a Galactic 
dust-to-gas mass ratio and size distribution. The IR SED also 
depends on the relative abundance between VSGs and BGs. 

We defer the discussion about future observations needed to 
disentangle the cloud structure and the dust size distribution to 
sect. 8. Also note that we cannot estimate the total dust extinction 
from the models, given the lack of observational constrains in 
the FIR. The Ay value may be higher than the value derived 
from our modeling if the total H2 column density including cold 
molecular gas is higher than that derived from the mid-IR H2 
rotational lines. If we assume that the molecular gas is diffuse 
and A^h = 2 x 10 20 cirr 2 , this corresponds to Ay =0.1. This is 
smaller than the average value derived from optical observations 
(Ay ~ 0.6 in the main shock region for the diffuse emission). 



7. Dust processing in the Stephan's Quintet shock? 

The Stephan's Quintet galaxy-wide collision is an extreme en- 
vironment where observations and modeling dust emission may 
provide insight into dust processing in shocks. The galaxy col- 
lision must have triggered shocks across the whole ISM. As 
discussed in Guillard et al. (2009), the shock velocity depends 
on the preshock gas density. Gas at preshock densities njj > 
0.2 cirr 3 has been shocked at velocities small enough (V s < 
200 km s -1 ) to cool, to keep most of its dust, and to become 
molecular within a few million years. To account for the H2 
emission, the molecular gas must be processed by low-velocity 
(5 - 20 km s -1 ) MHD shocks, repeatedly. Therefore, the origin 
and dynamical state of the SQ molecular gas is very different 
from that of the Galactic ISM. 

In the previous section, we assume that the dust properties 
in the SQ ridge are identical to those of the Galaxy, which obvi- 
ously is a simplifying assumption, which may be far from real- 
ity. In the shock region, various processes can affect dust grains 
(e.g. thermal sputtering of grains in the hot gas or destruction in 
shocks due to gas-grain and grain-grain interactions, see Jones, 
2004, and references therein) that can affect both the dust-to-gas 
mass ratio and the dust size distribution. So far, observational ev- 
idence of dust processing comes mainly from gas phase metals 
depletions (Sembach & Savage, 1996). To our knowledge, there 
is no direct observational evidence for changes in the dust size 
distribution that can be unambiguously associated with shock- 
processing. Stephan's Quintet is an outstanding target to look 
for such evidence on galactic scales. 
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Fig. 8. Comparison between the two models presented in this paper: diffuse (dashed blue line, same model as Fig. 7) and clumpy 
(black line, same model as Fig. 7) molecular gas. The overlaid points indicate the Spitzer fluxes extracted within the shock region 
(see bottom right inset). The gas lines are removed from these spectra, as well as the IR component of the incident ISRF (Guv = 1-4). 
For comparison, the clumpy model with stellar continuum is shown on the thin black dotted (same model as Fig. 7). 



Dust destruction processes depend on the type of shocks. 
Guillard et al. (2009) show that within the timescale of the 
galaxy collision (~ 5 x 10 6 yr), the destruction of grains smaller 
than ~ 0. 1 /im is complete in the hot, teneous gas (corresponding 
to pre-shock densities «h 5= xl0~ 2 cirT 3 ) that is shocked at high 
velocities (V s > 300 km s~'). For intermediate preshock densi- 
ties (« H ~ 0.2 - 0.01 cnr 3 , y s = 100 - 300 km s" 1 ), models 
predict significant dust destruction (10 - 50 % in mass) and pos- 
sibly the production of an excess of small grains by shattering 
(e.g. Jones & Tielens, 1994). Within the dense («h > 10 3 cm -3 ) 
molecular gas, low-velocity MHD shocks (5 - 20 km s 1 ), may 
only have little effect on dust (Guillet et al., 2007; Gusdorf et al., 
2008). 

The present SQ data provide some constraints on the dust 
size distribution. The detection of PAH emission from such a vi- 
olent and extreme environment as the SQ shock may be seen as 
a surprise. PAHs are predicted to be completely destroyed for 
V s > 125 km s _1 («h = 0.25 cirT 3 ) and their structure would 
be deeply affected for 50 < V s < 100 km s -1 (Micelotta et al., 
2010). Therefore, the PAH detection provides interesting con- 
straints on the density structure of the preshock gas. Since it is 
unlikely that PAHs reform efficiently from carbon atoms in the 
postshock gas, we conclude that (i) they were protected from 
high-velocity shocks in high preshock density regions (nn S 
0.3 cm -3 ), or/and ( ii) they are the product of the shattering of 
VSGs in the shock. 

We note that the 6.2 pm AIB is absent and that the 17/itn 
complex is prominent. This suggest that large PAHs, which emit 
more efficiently at larger wavelengths (e.g. Draine & Li, 2007), 
would be predominant in the shock. This may result from PAH 



processing in shocks, larger molecules being less fragile than 
smaller ones (Micelotta et al., 2010). This interpretation is sup- 
ported by Spitzer observations by Tappe et al. (2006) that show a 
prominent 17 pm emission from the supernova remnant N132D 
in the Large Magellanic Cloud. 



The 12/24 //m flux ratio is sensitive to the relative mass 
abundances between PAHs and VSGs. Fig. 9 shows the pre- 
dicted 12/24/im flux ratio as a function of the PAH / VSG dust 
mass ratio. For each PAH/VSG ratio, the dust emission spec- 
trum is calculated for an ISRF of Guv = 1.4 with the DUSTEM 
code. We derive the 12 and 24pm fluxes by integrating the 
model spectrum over the 12 pm IRAS and 24 pm MIPS filter 
bandpasses. The PAH/VSG mass ratio in the SQ shock is re- 
markably close to the value for the Galactic diffuse ISM. A de- 
viation would be expected if the dust had been processed by 
high-speed (> 100 km s _1 ) shocks (Jones et al., 1996). Like for 
PAHs, this may suggest that the postshock dust was lying in gas 
dense enough to have been protected from destruction by fast 
shocks, and thus was protected from the effect of fast shocks, 
as proposed by Guillard et al. (2009). However, as discussed in 
sect. 4.3.2, the PAH spectrum in the shock is significantly differ- 
ent from that observed in the diffuse Galactic ISM. The fact that 
the long wavelength bands (11.3 and 17yum) are brighter than 
the 6.2 and 7.7 pm features imply that a higher fraction of the 
PAH emission is emitted within the 12 pm IRAS band. If true, 
this correction would imply a lower PAH/VSGs abundance than 
in the Galactic diffuse ISM. 
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Influence of PAH/VSG ratio on 12/24/im flux ratio 
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Fig. 9. Predicted 12/24 //m flux ratio from the DUSTEM model 
as a function of the PAH / VSG mass ratio. The horizontal red 
and blue lines shows the observed 12/24 //m flux ratios for the 
SQ shock, and the Galactic diffuse ISM (IRAS observations). 
The vertical blue line is the Galactic PAH/VSG ratio adopted in 
the model (0.915). This value is chosen to fit both the 12/24 /mi 
ratio and the extinction curve. 

Far-IR SED imaging of the dust, possible with the PACS 12 
and SPIRE 13 instruments onboard Herschel, combined with 
modeling of the full IR spectral energy distribution, would pro- 
vide diagnostics for measuring the relative abundances of the 
different grain populations. 

In our modeling (sect. 6), we did not consider any thermal 
emission from collisionally heated dust in the hot plasma. This 
emission may arise after a fast shock wave has gone through 
tenuous, dusty gas, but it will last a very short period of time 
(~ 10 6 yr), producing a "flash" of FIR emission, because the 
dust cooling efficiency drops as the grain sputtering occurs in 
the hot (> 10 6 K) gas (Smith et al., 1996; Guillard et al., 2009). 
However, there are two reasons why this may not be a valid 
assumption. First, if there is a significant dust mass in grains 
larger than about 0.3 /mi, this dust may survive for > 5 x 10 6 yr. 
These grains may contribute to the FIR emission, as proposed by 
Xu et al. (2003) and further investigated by Natale et al. (2010). 
Second, we cannot exclude that some dust may be injected 
into the hot phase by the ablation of clouds (turbulent mixing), 
due to their dynamical interaction with the background plasma 
(Guillard et al., 2009). A continuous supply of dust from the 
warm to hot phase could balance destruction by sputtering. This 
possibility needs to be quantified, which is beyond the scope of 
this paper. 

8. Summary and concluding remarks 

In this paper, we present new Spitzer imaging and spectro- 
scopic observations that reveal PAH and VSG emission from the 
galaxy-scale shock in Stephan's Quintet (SQ). Here are our main 
observational results: 

- Faint dust emission is detected in the center of the H2-bright 
SQ shock structure, outside star-forming regions lying in the 

12 Photodetector Array Camera & Spectrometer 

13 Spectral and Photometric Imaging Receiver 



SQ halo. The 12/24 /mi flux density ratio in the SQ ridge is 
remarkably similar to that of the diffuse Galactic ISM. This 
suggests that the PAH to VSG abundance ratio is similar to 
that of the diffuse ISM of the Galaxy. 

- The global mid-IR SED is consistent with the expected dust 
emission from the amount of warm H2 detected by Spitzer 
(Nh — 2 x 10 20 cm" 2 ) for a UV radiation field intensity of 
Guv ~ 1 [Habing unit], which is consistent with UV obser- 
vations of the shock. 

- The PAH emission spectrum in the SQ shock is significantly 
different from that of the diffuse Galactic ISM. The 7.7, 1 1.3 
and 17 /mi aromatic bands are detected, but the 6.2 /mi band 
is absent. The 17 /mi complex is prominent, but the 16. 4 /mi 
is not detected. Interestingly, the 7.7/11.3/mi flux ratio in 
the SQ shock is a factor ~ 2 lower than that of the dif- 
fuse Galactic PAHs. These characteristics may suggest an 
enhanced fraction of neutral and large PAHs. 

Spitzer imaging and spectroscopy reveal powerful H2 emis- 
sion in the Stephan's Quintet X-ray giant shock that extends over 
the full area (* 35 x 15 kpc 2 ) of the ridge (Cluver et al., 2010). 
We expect some dust emission to come from molecular gas be- 
cause H2 forms on dust grains. In this paper we test this inter- 
pretation by modeling the IR emission from dust associated with 
the H2 gas present in the SQ shock structure, and by comparing 
the model calculations with Spitzer observations. 

- We model the emission from dust associated with diffuse or 
clumpy molecular gas, embedded within Hn gas and X-ray 
emitting plasma. The model SED is consistent with mid-IR 
Spitzer observations for both cases, for a Galactic dust-to-gas 
mass ratio and a Galactic dust size distribution. For diffuse 
gas, the best-fit column density is A^h = 1.8+0.5xl0 20 cirT 2 , 
which is close to the value derived from warm H2 observa- 
tions. For clumpy molecular clouds that are optically thick 
to UV radiation, we find that the H2 surface filling factor is 
/ ~ 0.03. So far, the present data and the degeneracy be- 
tween the dust size distribution and the cloud structure do 
not allow us to decide whether the molecular gas is diffuse 
or fragmented into clouds that are optically thick to UV pho- 
tons. 

- The presence of dust in the SQ shock shows that dust is able 
to survive in such a violent environment. We propose that 
at the time of the high-speed galaxy encounter, the dust that 
survived destruction was in pre-shock gas at densites larger 
than a few O.lcrrT 3 , which was not shocked at velocities 
larger than ~ 200 km s (see Guillard et al., 2009). Present 
data do not allow us to indentify a plausible impact of the 
shock on the dust size distribution, but the peculiar proper- 
ties of the PAH emission in the SQ ridge (summarized above) 
may be the result of PAH processing in shocks. 

Future Far-IR observations are needed to constrain both the 
structure of the H2 gas and the dust properties in the SQ shock. 
The PACS and SPIRE instruments on Herschel would provide 
SED of the Far-IR dust emission resolved on scales of 5 - 10 arc- 
sees, crucial for helping to distinguish how clumpy the dust dis- 
tribution is in the shock structure. The flux in the FIR is different 
for the two models (diffuse or clumpy, see Fig. 8). The spatial 
resolution of the PACS [60 - 210/mi] instrument on-board the 
Herschel Space Observatory at 70 pm is comparable to that of 
Spitzer at 24 pm. It will thus be possible to extract the FIR emis- 
sion from the shock with much more accuracy than with Spitzer. 
This is needed to possibly break some degeneracies of the model 
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and better constrain the key-parameters (e.g. dust size distribu- 
tion). 

The Stephan's Quintet galaxy-collision is a unique environ- 
ment to study dust survival in shocks, on galactic scales. The 
presence of dust in the SQ shock stresses the need to revisit the 
standard calculations of dust survival timescales in the ISM, by 
taking into account its multiphase structure. 
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Appendix A: Fitting PAH features 

The results of the PAHFit decompositions of the Spitzer IRS 
spectrum of the center of SQ ridge and, for comparison, of the 
ISO-CVF spectrum of the diffuse Galactic medium, are pre- 
sented in Fig. A.l, A. 2 and A.3. No extinction parameter is 
introduced in the fit. 
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PAHFit: IRS Spectrum of SQ Shock 
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Fig. A.l. Result of a PAHFit run on the 5-38 pm Spitzer IRS spectrum extracted over a ~ 18" x 15" area centered on the position 
in the center of the ridge (ON beam on Fig. 1). The green line shows the complete model. The fitted gas lines are in magenta. The 
gaussian profiles of the gas lines are used to remove their contributions and extract a "pure dust spectrum" (see Fig. A. 2). 




Fig. A.2. PAHFit decomposition of the Spitzer IRS dust spectrum (smoothed over 5 resolution elements, and from which gas lines 
have been removed) extracted in the center of the SQ ridge. Blue solid lines shows the lorentzian components of the PAH decom- 
position, and the thick gray line is the total (dust + starlight) continuum. The result of the fit is the green solid line. 
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Fig. A.3. PAHFit decomposition of the ISO-CVF spectrum of the diffuse Galactic light (Flagey et al., 2006), centered on the Galactic 
coordinates (26.8, +0.8). The gas lines have been removed from the spectrum. The lorentzian components of the decomposition of 
the PAH features are shown in blue. All components are diminished by the extinction, indicated by the dotted black line, with axis 
at right. The solid green line is the full fitted model, plotted on the observed flux intensities and uncertainties. 



